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This paper presents a combination between a fuzzy logic control (FLC) and a 
predictive direct power control for multifunctional grid connected photovoltaic 


(PV) system, to solve the oscillation problem in the DC link voltage of the in- 
verter caused by the fast irradiation changing. The whole system consists of a PV 
system which interface a DC-AC inverter, a FLC maximum power point tracking 
(MPPT) algorithm has been adopted to operate the DC-DC converter at the MPP. 
The predictive control strategy is applied to the DC-AC inverter with FLC in its 
voltage control loop to improve the power exchange between the grid and the PV 
system. Simulation results have been verified through MATLAB/Simulink soft- 
ware for the purpose of giving the effectiveness of the suggested control against 
existed controllers. 
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1. INTRODUCTION 

Recently, with the rapid depletion of traditional sources of energy, renewable energy source has be- 
come the mainstream of industrial applications [1]. Among them, solar energy is the most promising one and 
has usually been in the best location to reduce the fossil fuel dependence in the energy industry [2], [3]. Par- 
ticularly, photovoltaic (PV) power generation, is conducting to a fast growing integration of renewable energy 
sources as a distributed generation (DG) system into the utility grid. However, a large scale penetration of these 
sources poses serious problems, as it can threat the power system’s stability and reliability [4]-[6]. Meanwhile, 
the main challenge of theses systems is the extraction of the maximum power of the PV system through max- 
imum power point tracking (MPPT)s algorithms, as well as, the power quality injected to the grid, which is 
in direct relation with power converter and its control strategy [7], [8]. Therefore, the performance of theses 
systems is highly dependent on the choice of the power converter and the appropriate control technique [9]. 
The DC link capacitor voltage of the distributed DC-AC inverter should be controlled through a robust con- 
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trollers in order to enhance dynamics, which is prone to rapidly changing of climatic conditions as well as the 
power load variations. To cope with this problem, several control strategies have been reported in the literature 
[10]-[15]. Lu et al. [16] proposed an enhanced model predictive direct power control for a grid-connected 
converters where, two voltage vectors are applied over a control period and their duty to optimize the benefit 
of the PV owners and satisfy the limits of the voltage variations simultaneously. Krommydas and Alexan- 
dridis [17], a new approach of analysis is employed on the basis of nonlinear systems theory with purpose to 
develop an appropriate controllers with guaranteed stability. Researchers in [18], [19] propose an optimized 
direct power control and a predictive current control respectively in grid connected inverter combined with a 
PV system through PI controller to improve the power quality generation as well as the whole system stability 
under rapidly environmental conditions variation. Wu et al. [20] proposes a dc voltage droop control in the 
grid-connected photovoltaic system with a view to make the capacitance of the medium time scale contribute in 
the grid frequency response without any supplementary equipment. Zhou et al. [21] proposes a robust DC-link 
voltage control based on the linear active disturbance rejection control scheme to enhance the tolerance of PV 
grid-connected inverter to disturbances. In this work, fuzzy logic control as a robustness controller in predictive 
direct power control is suggested to enhance dynamics and provide a good quality of power to the grid in the 
presence of non linear load under all dynamic conditions. The purpose of this work is to propose a suitable 
direct power control (DPC) scheme combined with a robust fuzzy logic control (FLC) controller for the multi- 
functional grid-connected PV system based on a predictive approach. The principle of this approach is based 
on the selection of the optimal control vector to be applied during the sampling period. Where, the selection 
is done by optimizing a performance function. Moreover, in the dc link voltage control a FLC controller has 
been proposed to enhance the problem of PI controller under external disturbances (fast changing irradiance, 
non linear load variation). 


2. SYSTEM DESCRIPTION 

The schematic diagram of the multifunctional grid connected PV system with its control algorithm 
scheme is depicted in Figure[I| The PV system (part 1) consists of a PV panels, a DC-DC boost converter which 
is interfaced the DC link bus of the multifunctional inverter. A FLC MPPT algorithm is exploited to maximize 
power, and to provide pulses through PWM for controlling the DC-DC boost converter. The multifunctional 
inverter (part 2) is connected to the grid at the point of common coupling (PCC), which is controlled through a 
predictive direct power control via a FLC regulator in its voltage loop as depicted in the scheme of the control 
algorithm. 


2.1. PV system 

The PV array is designed in this paper for a 3600 W peak power tracking capacity. A solar PV module 
kyocera KC 200 W adopted in this work, the PV generator is composed of three strings, each string contains 
6 panels. In order to interface the dc link voltage bus, the output voltage of the DC-DC converter must be 
very high. Hence, a DC-DC boost converter is necessary to offer high voltage gain as well as to reach the 
maximum power point via MPPT algorithm, the duty cycle D is the output of the MPPT algorithm which will 
be investigated to produce control signal for the switching metal oxide semiconductor field effect transistor 
(MOSFET). The PV system control scheme block is depicted in Figure [I] highlighted (part 1). The DC link 
voltage bus is given as (1): 


Vow 
Vde = 2 


=i p (1) 


2.1.1. Maximum power point tracking fuzzy logic control 

The investigation of fuzzy logic control has become more attractive during the last decade since it can 
deal with imprecise inputs, doesn’t need any accurate mathematical model. Therefore, the application of FLC 
as a MPPT algorithm provides a good performance, where, the MPP tracking doesn’t require any parameter 
information and consists of a step-wise adaptive search which leads to a rapid convergence [22]. This algorithm 
purposes to reach the MPP from the PV array whatever climatic conditions. The controller consists of three 
stages namely fuzzification, rules inference and defuzzification as presented in Figure[2] 
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Figure 1. The schematic predictive DPC control of the whole system 


Crisp input Crisp output 


Fuzzification Rule Inference Defuzzification 


Figure 2. Structure of FLC 


In this article, Ap(k) and Av(k) are the input variables at each sampling time k of the proposed algo- 
rithm, the output variable is D which is the changing of the duty cycle. The variable Ap(k) and Av(k) will be 
given as: 


Ap(k) = p(k) — p(k — 1) Av(k) = v(k) — v(k — 1) 


where p(k) is the PV array power and v(k) is the PV array voltage. Thus, at the MPP Ap(k) and 
Av(k) equal zero. Over fuzzification, the variables Ap(k) and Av(k) are transferred into linguistic variables 
through a membership functions. The Ap(k) and Av(k) universe of discourse are adopted to five fuzzy sets, 
including negative big (NB), negative small (NS), zero (ZE), positive small (PS), and positive big (PB) as 
depicted in Figure B] Then, the output variable is developed based on rule base presented in Table 1. 
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Ap(k), Av(k) 
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Figure 3. Membership functions of: Ap(k), Av(k) and A(k + 1) 


Table 1. Rule base table of FLC MPPT 
Av(k) 


Ap(k) ZE ZE ZE ZE ZE ZE 


Mamdani fuzzy inference process is investigated with a fuzzy max-min operation in order to provide 
the fuzzy algorithm MPPT output. The obtained crisp is calculated in the defuzzification step via the center of 
gravity method. It computes the center of gravity of the final fuzzy crisp. The output value of the algorithm 
AD(k) is provided as (2): 

AD(k) = Sa MAD: (2) 
dino Wi 


where W; is the weighting factor, n is the maximum number of effective rules and D; is the value corresponding 
to the membership function of D. Finally, the duty cycle will be expressed as D(k + 1) = D(k) + AD(k) 


2.2. Multifunctional inverter 
2.2.1. Predictive direct power control 

From Figure |4| the influence of filter and source resistors is omitted, the system equations are ex- 
pressed as (3)-(5): 


dis 


e — Va = Larry (3) 
diş 

Uf — Vs = Lf i (4) 

istip =i (5) 


The purpose of the multifunctional system is for obtaining a sinusoidal current of the utility is, using 
the relationship [5] Assume that 7; can be expressed as (6): 


=s + tn = ts + tf (6) 


where i; is the fundamental current and 4; p is the harmonics currents. Then, by manipulation in (6) can be 
exhibited as (7): 


11, f — ts = — (ilh = if) => Ais = —Aiş (7) 
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Figure 4. Schematic circuit of the multifunctional inverter 


Based on the approximation Ai = di for a small current variation; (7) will be expressed as (8): 


dis = —dis (8) 
By subtracting (4) from (3) and substituting (5) in the result equation, we will obtain: 


di 
up —e= (Ly + Ls) = (9) 


The obtained equation depicting the multifunctional inverter is presented as (10): 
uf — e= L— (10) 
where L = Ly + Ls 


In the stationary a-( reference frame and for a balanced three-phase system, the dynamics of the currents 
drawn by a multifunctional inverter [23], is given as (11): 


ji 


dija — L (ey — Va) 
dt T © \a— Ya (11) 
z (es — ve) 


From this equation, a discretization of the first order, over a switching period T,, the components of the current 
vector at the end of the switching period are given by (12): 


Aia = talk + 1) — talk) = 7 (ea(k) — valk) 
(12) 
Aig = ig(k + 1) — ig(k) = + (ea(k) — va(k)) 


Furthermore, the instantaneous active and reactive powers are expressed by (13): 


P |_| €a eB to 
Pa ENI a» 


The sampling period T, is supposed to be small in comparison with the period of the power-source 
voltage, components of eag are supposed constant during the switching period ea(k) = ea(k + 1). Conse- 
quently, the variation of active power Ap and reactive power Aq can be given as (14): 


Ap = €g(k) - Aia + eg(k) - Aig 
(14) 


Aq = eg(k) - Aia + €a(k) + Aig 


Just seven variations can be realized on the active and reactive powers, by applying each of the seven 
control vectors over a switching period. These variations, named Ap;(k) and Aq;(k) are given by (15): 


Api = Ẹ leas (RII? — % lleas(k)I| + [cos()vas + sin(A)u gi] T 
Agi = F lleag(K)|| - [cos(0)vgi + sin(4) vai] 
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In normalized quantities these variations are expressed by (16): 


Ne Api = lleas ll =. s =. 
Api = Ta Jjeasliloaall ~ Ivas [cos(0)Uai + sin(0)Ug:] 
(16) 
Aq = r Aq = cos(0)Ugi + sin(@)ta; 1=0.1....,6 
T leasllllvasll 
The desired variations, Ap;(k)* and Aq;(k)* , can be written in normalized quantities as (17): 
Api = AP o 
F lleasllllvasll 
(17) 


aS Aq; 
Ti = Flleaallllvaall 
2.2.2. Prediction of instantaneous power references 
The P-DPC needs the prediction of the instantaneous active and reactive power references at the next 
sampling time p*(k + 1) and q*(k + 1). Reference reactive power is frequently imposed zero for operation 
with a unitary power factor. The estimation of the reference of this power at the next sampling time (k + 1) is 
given as (18): 


q (k +1) = g” (k) (18) 


The active power command p*(k) is obtained from the FLC dc link voltage controller. If tracking error of 
dc link voltage is supposed constant during two successive sampling periods, the instantaneous active power 
command p*(k + 1) can be estimated through a linear extrapolation as (19): 


p(k +1) = 2p* (k) — p*(k — 1) (19) 


2.2.3. Principle for selecting the optimal control vector 

This predictive approach is based on the application of a single control vector at each switching period, 
which is similar to the DPC based on a predefined switching table, except that the control vector applied at this 
time is optimal. It corresponds to the most adequate variation among the seven possible variations described by 
(17). Once the desired variation [Api(k), Aqi(k)] is calculated and normalized, the optimal control vector to 
be selected over the switching period [kT;, (k + 1)T] corresponds to the variation point [Api ; Aq] closest 
to the desired variation point. Then, the performance function to be optimized consists in searching for the 
smallest distance between the point [Ap Agi" and the points [Ap;(k), Aq;(k)]. It is formulated by the 
following objective (20) [24]: i 


o 12 N 2 
f= min (fe Api! [Aq aal’) i=1,....6. 20) 


2.2.4. Fuzzy logic controller for DC link control 
To elaborate the fuzzy logic controller, the input signals for this controller is the error of dc link 
voltage and its variation. While, its output signal is the variation of the maximum value of the current Asmax 


as presented in Figure[5] [25]. 
Decision Macking 


Alsmax 


Figure 5. Schematics of dc link voltage regulation through FLC controller 
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The fuzzy controller is developed as: seven fuzzy sets for each input and output: negative big (NB), 
negative medium (NM), negative small (NS), zero (ZE), positive small (PS), positive medium (PM), positive 
big (PB). Mamdani method is employed with a fuzzy max-min operation. The membership functions of error 
(€) and change in error (de) of dc link voltage are shown in Figure 6. 


€, de 


NB NM NS PS PM PB 


-1 -0.5 -0.25 0 0.25 0.5 1 


Figure 6. Normalized triangular membership functions 


2.2.5. Design of rules table 

The fuzzy control rules which relates the input variables to the output are presented in Table 2. The 
deffuzification process produces the single crisp value of output (Asmax) from the aggregated fuzzy set which 
includes a range of output values. The commonly used centroid method is employed to transform the fuzzy 
subset of (Asmax) to numerical value. It calculates the center of gravity from the obtained output fuzzy set, 
and provides the final result. 


Table 2. Rule base table of FLC controller 


NB NM NS ZE PS PM PB 

NB NB NB NB NB NM NS ZE 

NM NB NB NB NM NS ZE PP 

NS NB NB NM NS ZE PS PM 

de ZE NB NM NS ZE PS PM PB 
PS NM NS ZE PS PM PB PB 

PM NS ZE PS PM PB PB PB 

PB ZE PS PM PB PB PB PB 


3. RESULTS AND DISCUSSION 

The proposed predictive direct power control combined with FLC controllers was verified in simula- 
tion using MATLAB/Simulink. The topology of the simulated system is depicted in Figure[I]and the simulation 
parameters of the entire system are given in the Table 3. In order to illustrate the simulation steady and dynamic 
state to assess the effect of fast irradiation changing, a consecutive irradiation ramp evolution was adopted. This 
irradiation change starts from 0 W/m?2, stops at 350 W/m2, waits at this level for 0.25 s, and then continue to 
increase again to 600 W/m? and 1000 W/m? respectively with a constant slope and stay at their levels during 
1.5 s as depicted in Figure[7] The environmental temperature is assumed constant (25 °C) in all simulation test. 


Table 3. PV array parameters 


Parameters Values 

Us, fs 80 V, 50 Hz 

Ls, Rs 0.1 mH, 0.1 Q 

Li, Ri 10 mH, 20 Q 

Lf, Rf,C 1mHĦ,0.01 Q, 1100 u F 
Le, Re 0.55 mH, 0.01 Q 
PVarray Npp=3, Nss=6 


Figure[8|depicts waveforms of the source current, the load current and the filter current and their zoom 
at each ramp irradiation changing. Figures [9|and[10|show the power flows as well as the duty cycle of the dc 
dc boost converter. When the control is not activated with a G=0 W/m? between 0 s and 0.1 s, non linear 
load causes a serious problem of current harmonic and reactive current at the pcc exhibited by the waveform 


An enhanced control scheme for multifunctional grid connected PV system ... (Boualem Boukezata) 


860 og ISSN: 2302-9285 


of the source current and the reactive power supplied by the grid. When the predictive direct power control 
is switched on at t=0.1 s, the source current becomes sinusoidal improved by the injected filter current with a 
unity power factor. The source current THD decreased to 1.66% as shown in Figure [1] which is acceptable 
with IEFE Standard 519-1992. 


Irradiance 


0 0.2 0.4 0.6 0.8 1 
t(s) 


Figure 7. Profile of irradiance 
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Figure 8. Waveforms of the whole system currents 
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Figure 9. The entire system power behavior 
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Figure 10. Duty cycle behavior 
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Figure 11. Grid current FFT 


Furthermore, in the multifunctional PV system mode at t=0.3 s, an increase in the irradiation at 0.3 
s gives rise a PV power generation, the PV system and the grid provides the required power demanded by 
the load. Thus, for each step irradiation changing at 0.6 s and 0.85 s respectively, the FLC MPPT algorithm 
provides a best tracking of PV power with a minimum oscillation; where there is no power loss in dynamic 
and steady state cases. For each PV power increase the active power fedded by the grid decreased until the 
multifunctional PV system supports all non linear load power demanded. The source current stays sinusoidal 
in phase with the voltage source whatever the irradiance variation, the THD level is reduced to 1.25%, as given 


in Figure 


fi 
09 0.002 0.904 0.006 0.008 091 a.012 0.014 aoe ome 


Fundamental (50Hz) = 19.43, THD= 1.25% 
T T T 


oe oe ee ee .E. B. „m 


s00 000 


Frequency (Hz) 


Figure 12. Grid current FFT 
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The dc link voltage follows its desired value 300 V at 0.1 s with a minimal overshoot as well as a rise 
time of 0.0015 s with FLC controller, and then for each step changing irradiation, the FLC controller proves 
more robustness and stability under fast irradiation changing; there is no overshoot just a little increase of DC 
link voltage then the controller can maintains it at its desired level, as shown in Figure [13] However; through 
PI controller; we show that the DC link voltage reaches its desired value before 0.165 s to reach 300 V with 
an overshoot of 28v. The same remark when the irradiation occurs at 0.3 s; 0.6 s and 0.8 s the increase in 
irradiation causes an increase in the dc link voltage with overshoot of 12.05 V then it reaches its level during 
0.1 s. 


Figure 13. DC-link voltage behavior under ramp disturbances changing in solar irradiation 


3.1. Comparative study between previously developped algorithms 

Comparative study between previously developed algorithms is simulated with the same parameters 
of the proposed algorithm and is given in Table 4 under different irradiance conditions, specifically G=0 W/m? 
and G=1000 W/m?. The THD levels are measured for each algorithm in both scenarios. In the case of Classical 
DPC, THD levels are recorded at 2.32% under no irradiance (G=0 W/m?) and 1.43% under irradiance (G=1000 
W/m?). The Improved DPC algorithm shows notable improvement, with THD levels reduced to 1.46% and 
1.3% under G=0 W/m? and G=1000 W/m?, respectively. The proposed predictive DPC algorithm demonstrates 
competitive performance, yielding THD levels of 1.66% and 1.25% under the respective irradiance conditions. 


Table 4. Comparative study of differents developped controls 


Algorithms THD level (%) THD level (%) 
G=0 G=1000 W/m? 
Classical DPC THD=2.32 THD=1.43 
Improved DPC THD=1.46 THD=1.3 
Predictive DPC THD=1.66 THD=1.25 


4. CONCLUSION 

Robust control algorithm in dc link voltage is very required in multifunctional grid connected PV 
system, which is directly affect the dynamic of the system. Therefore, a FLC controller with predictive direct 
power control was proposed in this system. A comparative study between direct power control algorithm with 
classical and new switching table and proposed predictive direct power algorithm have been verified in term of 
THD level. Simulation results confirm the superiority of the proposed controller against previously developed 
controllers. 
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